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1. Introduction 

The fascinating progress in the study of graphene [1, 2] in many ways is based on the 
usage of various experimental methods that provide new information on the intriguing 
electronic properties of this material. An impressive example is provided by the 
measurements of ac response which includes reflectance, transmission and extracted 
from them conductivity of monolayer graphene [3-5], bilayer graphene [6], epitaxial 
graphene [7], and graphite [8] reported in 2008. These measurements have made possible 
a quantitative comparison of the value of the universal ac conductivity background 
extracted from the experimental data and the theoretical result [9, 10] predicted earlier 
which at zero temperature can be written in a particularly simple analytical form 

Here Q{> 0) is the photon energy, is the chemical potential tunable by the gate voltage 
and r is the impurity scattering rate which is set to zero in the second, interband, 
term. It is remarkable that the interband conductivity is universal, independent of the 
band structure parameters. The dependence of the absorption edge on concentration in 
semiconductors is called as the Burstein-Moss effect. It allows one to measure directly 
the value of the chemical potential of the system. For finite temperature versions of 
Eq. ([T]) see, for example, Refs. [10-12]. 

The theoretical value vre^/ {2h) of the universal ac conductivity of monolayer 
graphene is in good agreement with experiments in the infrared range of wavelengths 
A ~ 1500 - 2500 nm [5] and even in the visible with A ~ 400 - 700 nm [4]. The last 
result seems to be well beyond the validity of the Dirac approximation used in [10, 11], 
but as shown in [13] the corrections to the Dirac-cone approximation are a few percent 
only. 

Even earlier, similar measurements were made in an applied magnetic field on 
epitaxial graphene [14, 15], monolayer graphene [16, 17], bilayer graphene [18] and 
highly oriented pyrolytic graphite [19]. The ac conductivity was extracted from the 
reflectance in this last paper only, while in [14, 16, 18] the relative transmission is 
considered. Knowledge of relative transmission allows one to investigate the magnetic 
field dependence of the energy levels and find a band velocity, although as yet does not 
allow one to study the redistribution of the optical spectral weight between different 
energy regions as done in [8]. On the theoretical side, the magneto-optical response of 
graphene has been studied in Refs. [20-23], and including the role of interactions and 
collective excitations, in Refs. [24,25]. 

Besides the agreement with early theoretical predictions [9-12, 20, 21], several 
observations [5] indicate the relevance of many-body interactions to the electromagnetic 
response of graphene. Those include unusual large broadening around 2|/i| which cannot 
be explained by thermal smearing of the curve described by ([1]). Also the conductivity 
below the threshold of 2|/i| does not vanish due to Pauli blocking as described by Eq. ([T]), 
but shows a significant residual conductivity as strong as 0.37re^/ (2/i) in a wide range 
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of energies at finite doping. 

An attempt to explain these results taking into account effects of disorder and 
phonons (both optical and acoustic) was made in Refs. [26,27]. Earlier the influence of 
optical phonons on the ac conductivity was studied in [28] and recently corrections due 
to electron-electron interaction in undoped (/i = 0) graphene were discussed [29-32]. 

Thus one may conclude that the investigation of the ac response of graphene is 
already reaching the same maturity as, for example, for high-T^ superconductors (HTSC) 
[see Ref. [33], Sees. IV B and C, for a review], where the connection between conductivity 
and quasiparticle self-energy is exploited to extract information about the bosonic modes 
and their interaction with carriers. The case of graphene is, however, more complicated 
than that of HTSC, because one has to take into account the presence of interband 
transitions. 

In this paper we present some of our latest results on the ac conductivity 
accentuating the effect of various parameters such as temperature, impurity scattering 
rate, doping, magnetic field and of interactions such as the electron-phonon interaction 
on the value of the universal background. 



2. Model and general expressions for conductivity 

The frequency-dependent electrical conductivity tensor aa/si^) is calculated using the 
Kubo formula 

= — K^,{n + ^0) = ^ + , (2) 

where the retarded current-current correlation function is given by 

/oo 
rfte^^*nf^(t), u^pit) = -^eit)^r ipiUt), MO)]) , (3) 
■oo 

V is the volume (area) of the system, p = exp{—j3Ho)/Z is the density matrix with the 
Hamiltonian Hq in the grand canonical ensemble, /5 = 1/T is the inverse temperature, 
Z = Tr exp(— /3i7o) is the partition function, and Jq are the total paramagnetic current 
operators. The diamagnetic or stress tensor (tq,^) in the Kubo formula ([2]) is a 
thermal average of the diamagnetic part which is nonzero only when a tight-binding 
model is considered. Accordingly, the restricted optical conductivity sum rule (see e.g. 
Refs. [34, 35] for a review) in graphene has a form [36] 



oo 

.2^2 r ^2 



i I dnRea^^in) = M = J^^ [^^(e(k)) - n^(-6(k))] e{k), (4) 

— oo 

where e(k) is full band dispersion of graphene [37], np^u) = [exp((co' — fi)/T) + 
is the Fermi distribution function and the integration is done over the two-dimensional 
Brillouin zone (BZ). At large ~ t, where t is the nearest neighbor hopping parameter, 
the expression for a^xi^) has to be derived from a tight-binding model [36]. As studied 
in Ref. [27], for such large energies it is necessary to include in the analysis a hopping 
parameter for second-nearest neighbors with a value of the order of O.lt. 
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Still, one may consider the partial sum rules in the low-energy regime [38] based 
on the linearized continuum approximation for the tight-binding Hamiltonian. A big 
advantage of this approximation is that it allows one to obtain many transparent 
analytical results that elucidate the physics of graphene. Even more such insight 
comes when the low-energy description of graphene is mapped onto 2 -|- 1-dimensional 
electrodynamics [39] (see Ref. [40] for a review). In many cases, the effective low- 
energy QED2+1 description provides a very good starting point for the investigation 
of various transport properties making it possible to profit from the analogies with 
particle physics and to use its powerful field theoretical methods. In this case, the 
low-energy quasiparticle excitations are described in terms of a four- component Dirac 
spinor \E'^ = [ipK+As, i^K+Bs, i^K-Bs, i^K-As) ■ The spinor (with a given spin index s = ±) 
combines the Bloch states on the two different sublattices {A and B) of the hexagonal 
graphene lattice and with momenta near the two inequivalent Dirac points {Kj^ and 
K_) of the two-dimensional Brillouin zone. The low-energy Lagrangian density for 
quasiparticles can be written in a relativistic-like form as 

£ = ^ ^^{t, r) [^7° {hdt - tfi^) + thvpl^D, + ihvpl^Dy - A] ^,(t, r), (5) 

<T=± 

where '^^ = \I'j^7° is the Dirac conjugated spinor, the Fermi velocity vp = Stacc / i'^^) ~ 
c/300 plays the role of the speed of light c. The orbital effect of a perpendicular magnetic 
field B = V X A is included via the covariant derivative Da = da + {ie/hc)Aa, where 
a = x,y and the Zeeman energy is included via splitting fi^- = — o-fisB of the 
chemical potential fi with fiB = e/i/ (2mc) being the Bohr magneton. The 4x4 matrices 
7^^ furnish a reducible representation of the Dirac algebra. Here, we use the following 
representation: 

7° = fi (g) To, 7" = -if2 (g) Ta, (6) 

where the Pauli matrices and (as well as the 2x2 unit matrices Tq and tq) act on 
the valley {K^, and the sublattice {A, B) indices, respectively. This representation 
is derived from the tight-binding model for graphene and to make the treatment more 
general we also include a mass (gap) term with A (see e.g. review [40]). Accordingly, 
the Kubo formula ([2]) can now be calculated using the electric current density operator 
jait,r) that follows from the Lagrangian ([5]): 

Jait) = J d^rjait, r), jait, r) = j^^j^ = -evF J2 ^-(^' r)7°*a(i, r)- (7) 

All results in this paper are obtained neglecting the vertex corrections. In this case the 
calculation of the conductivity reduces to the evaluation of the bubble diagram. The 
corresponding polarization operator can be written in the form 

00 

- ' dujnF{uj) I ——- (8) 



V 2m J J (27r) 

—00 

X tr [7"(G^(cu, k) - G^{uj, \i))-i\G^{uj - 1], k) + G^(cu + fi, k))] 
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where the tr also includes the summation over spin index and G^^^\lj, k) are retarded 
(advanced) fermion Green's functions (GF). The approximation of Eq. ([8]) does not 
correctly reproduce the Boltzmann dc conductivity because of the omission of vertex 
corrections. The current vertex corrections for electron-phonon scattering in graphene 
within the Kubo formalism for the dc conductivity were considered in Ref. [41]. For 
the Coulomb interaction and ac conductivity the vertex corrections were considered in 
Ref. [32]. In the absence of magnetic field these GF are given by (we also set h = 1, so 
that fl can be regarded both as the frequency and energy of the photon) 

G->,K- -f'^^::'^'7^^^ - (9) 

In ([9]) we included the retarded (advanced) electronic self-energy T,^^'^^{uj) which 
represents the contribution due to impurities, phonons or electron-electron interactions 
and which in general also depends on the wave vector k. However, because we study self- 
energy corrections due to optical phonons in what follows we omit this dependence. We 
write all expressions in terms of the retarded self-energy S^(co') related to the advanced 
by a complex conjugation T,^{uj) = [S^(u;)]*, so that the subscript "R" will be dropped. 

The calculation of the conductivity follows closely that described in Appendix B of 
[21], so we present a final result for the case of zero magnetic field 

oo 

o2 



e f 

—oo 

(cu-S(cj))(cj'-S(cu')) - A2 



[s(cj') - j:{uj) -n][2uj + n- e{uj) - e(cu')] 

(o;-S(o;))(c^'-E*(a;')) - A2 



(10) 

[s*(.') - s(.) -n][2u + n- E(^) - j:*{u')]J " " ^(^^^'l + ^ ^'^ 

where u' = u + Q. The expression (ITOl) can be explicitly calculated as was done in 
Ref. [21] for ReS(cj) = and in Ref. [27] for A = including both ReS(cj) and ImE(cj). 
These final expressions are somewhat cumbersome as compared to the representation 
given in (fTOl) . 

3. AC conductivity in the absence of an external magnetic field and 
electron-phonon interaction 

As mentioned in Sec. HI for the correct analysis of the ac conductivity with phonons it 
is crucial to include both real and imaginary part of the self-energy S(c<j). However, to 
gain some first insight into Eq. (ITOl) it is useful to consider first the case ReS(co') = 0. 
Considering the imaginary part of the self-energy, r(ci;) = — ImS(u;) to be small, one 
can obtain the following representation [10] for the conductivity 

oo 

Rea..(fi,T) =^/ciu; ^ (11) 
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X 



2T(uj) 2T(uj) 



{\UJ\ + \UJ'\){UJ^+UJ'^) 



where the first term in square brackets of Eq. ffTTj) describes the intraband transitions 
and the second term describes the interband transitions. In the clean limit, F — >■ 0, both 
terms in the square brackets reduce to delta functions and we obtain a simple analytical 
expression, 



(12) 



Reaxx(^i) = — ;— t^t; In 2 cosh — - H — np \ — np { — ] 

^ ^ h n^ + 4T^ V 2TJ 2h [ \ 2 J \2 J 

where as in Eq. ([1]) the impurity scattering rate is only kept in the intraband (Drude) 
term for a regularization. In the limit F the real part of the intraband conductivity 
can be written Dirac delta-function 

Rea-*™(n) = -^In [2 cosh (A)] 5^. (13) 

The analysis of experimental data [5, 7] created a need for various generalizations of 
Eq. ([1]) to finite temperature, finite impurity scattering rate in the interband term and 
a nonzero gap A. Postponing the analysis of formulae with nonzero gap to the end of 
the section, we rewrite Eq. fll2p in a slightly different form (see also Refs. [11,12] and a 
review [42]), 

Re a-*™(fi) = ^ 7^^^^ In (2 cosh ^) , (14) 



hVL'^ + AT^ V 2T) 



sinh(^]/2T) 



2h cosh(/i/T) + cosh(fi/2T) ' 
For a finite T and small VL the interband term in Eq. f[T^ is linear in VL with decreasing 
slope as /i is increased. For i7 ^ and Q/T ^ 1 we recover the universal background 
value. So finite temperature effects do not deplete this background. These features 
are shown in Fig. [T^ where we show numerical results for the interband conductivity 
at two temperatures T = lOK and T = 50K for the neutrality point, /i = 0, and for 
/i = 40meV. We notice the depletion below 7re'^/2h of the low frequency region with 
increasing temperature and particularly with increasing value of the chemical potential. 
At large frequency Q the background remains unaltered. In all cases, however, the 
missing spectral weight is transferred to the intraband Drude term not shown in the 
figure. To establish this we show that a sum rule applies to the spectral weight under the 
real part of the conductivity. Integrating Rea^xi^) up to A chosen to be A ^ T 

we find 

^ Tre^ 

Rea^^(fi)(in = Wintra + Winter = (15) 

where the intraband and interband spectral weights are 

27re2 / UN vre^ /cosh A/2T + cosh(/i/T)\ 

W^ntra = lu (2 COsh -) , Winter = lu ( 1 + eOsh(^/T) ) ' ^'^^ 
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Figure 1: (a) The interband conductivity Kea^^'^^ (fl) in units 7[e'^/2h (universal 
background): for fi = OmeV the long dashed (red) line at T = 10 K, the solid (blue) line 
at T = 50 K, and for /i = 40 meV dash-dotted (black) line at T = 10 K, short dashed 
(green) line at T = 50 K. There are no interactions and the impurity scattering is set to 
zero, (b) The conductivity Kea^xi^) in units 7re^/2/i for /i = 40meV. The long dashed 
(red) line at T = 10 K, the solid (blue) line at T = 50 K both are for F = 2 meV. The 
dash-dotted (black) line at T = 10 K, short dashed (green) line at T = 50K are for 
F = 10 meV. 



For A ^ \fi\ T we arrive at the last equality in Eq. f|T5|) . This value is exactly equal 
to its zero temperature and zero /x value, when Reaxx{^) is simply equal to ne^ /2h for 
all Vt. Our results show that the interband contribution to the dc conductivity is zero for 
any finite temperature. However, there is also an intraband Drude contribution which 
in the clean limit gives an infinite contribution to the conductivity [see Eq. ( fT3l) ]. This 
agrees with the physical expectation that cr^c should be infinite at any finite T even 
at the neutrality point /i = 0. There are always thermally excited electrons and holes 
present at finite T and these do not scatter in the clean limit, so the conductivity will 
be infinite and not zero as claimed by Peres and Stauber [43]. 

In Fig. [lb we show results for the case of increasing impurity scattering F. In 
both cases shown the chemical potential is taken to be 40meV and two values of the 
temperature are used T = 10 K and T = 100 K. The long-dashed (red) curve and solid 
(blue) curve are for F = 2meV, while the dash-dotted (black) curve and short dashed 
(green) are for F = lOmeV. What is important to note is that for Vt > 200 meV all 
curves have risen to the universal background value, equal to one in our units. The 
width of the absorption edge at Q = 2fi = 80meV depends both on T and F. For 
example, for the same temperature larger F corresponds to a wider edge. 
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3.1. AC conductivity in the presence of a gap 



Next we consider the possibility of the opening of a gap A. The most general expressions 
for the complex conductivity which include a finite A were derived in Ref. [23] [Eqs. 
(13) and (14)] then rederived in Ref. [7] and used to analyze conductivity data in 
epitaxial graphene (see also Ref. [44] for an alternative scenario which preserves the 
linear dispersion law of carriers in graphene). In this case the clean limit generalization 
of Eq. dTl) is 



Recr! 



intra 

XX 



inter 



27Te? 
~h~ 



vre 



(5(fi) 



duj{—n'p{uj)) 
sinh(fi/2T) 



A2)^(CU2- A2) 



+ 4A2 



(17) 



e{Q^ - 4A' 



ef^xx V""/ 2/i cosh(/i/T) + cosh(^]/2T) 
In the T = limit the last expression reduces to a result given in Ref. [10] 



(yxx{^) 



A' 



+ 



9{Q-2msix{\fi\,A)). (18) 



h ' 2h 

The first term in Eq. f[T5]) gives the modified intraband Drude contribution and the 
second is the modified interband background. It is clear that for Q ^ |/u|,A this 
background again reduces to its universal value. 



4. Optical conductivity in the presence of electron-phonon interaction 

One may notice that the expression for the GF does not contain the chemical 
potential. This is because during the derivation of Eq. ([8]) the replacement oj ^ u — jj, 
moves /i from the GF to the Fermi function nF{oj), so that in the GF ([9]) the value 
uj = fi corresponds to the Fermi level. In the usual many-body calculations /i is kept 
in the GF. This implies that when the corresponding self-energy is inserted in the GF 
its dependence on energy has to be accordingly shifted u ^ u — fi. For example, 
for the simplified single-parameter model proposed in [45] that captures the qualitative 
features of the e-ph interaction in graphene, the imaginary part of the self-energy reads 

ImS(co', /i) = —G\uj — sgn(co' — ij,)ujo\ 9[{uj — — u^]. (19) 

Here ujq ~ 0.2 eV is the energy of the in-plane optical phonon mode and G is the 
dimensionless parameter determined by matching the model self-energy with the full ah 
initio result [45]. The corresponding real part of the self-energy can be obtained from 
Eq. ( |T9l) using Kramers-Kronig relations or directly calculating the self-energy for Dirac 
quasiparticles coupled to a single Einstein phonon [28,46,47] 

w'^{iJJ-^l+iJJo) 



ReS(c<j, ji) 





ujq In 




— uj\n 








uoq In 


{ul-^lf-uJl 





/i > 0, 
/i < 0. 



(20) 



{u)-uJo)'^{uj-^i+u}o) 

7eV is the energy cutoff which preserves the number of states. 



Here W = vW^ 

with t being the nearest neighbor hopping parameter. By definition, the carrier effective 
mass renormalization Ae// due the electron-phonon interaction is given by 



ReS(co', /i) = —Xeff{id — /i) + ReS^ 



[UJ 



(21) 
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where from fl2UI) we find that 



A.„.!^lf,„^_l + M, (22) 



IT \ l/i + CUol ^^0. 

The role of Ae// will be clarified in what follows. Choosing the value of G = 0.1 as done 
in Fig. [2^ below, for /i = one obtains Ag// ~ 0.16. This is of order of the theoretical 
estimate [45] and smaller than the value found experimentally [48,49]. 

So far we did not distinguish the bare chemical potential fi of noninteracting system 
from the chemical potential jl of the interacting system. The value of jj, is set by the 
doping p = sgn{fj,) fi"^ / irh'^Vp [here h is restored] which is controlled by the gate voltage. 
In the interacting system the dispersion law of the quasiparticle excitations is determined 
by the pole of the GF (jH]). For A = and neglecting the imaginary part of the self-energy 
one gets 

u — ReT,{u!, fi) = ±VF\ii\, (23) 

or counting the energy E of the excitations from the true Fermi surface of the interacting 
system, u ^ E + p, we obtain from Eq. (!23|) that 

E + fl-ReJ:{E + fi,fi) = ±VF\k\, (24) 

where in the vicinity of Fermi surface ReS(i? + fl,fi) = —XejfE + ReS(/i, fl) as follows 
from Eq. fl2T|) . Hence, we arrive at 



E{1 + Xeff) + fl-ReJ:{fl,jl) = ±VF\k\ (25) 

which allows us to identify the quantity 

fi — ReT,{fl, fl) = fi (26) 

as the bare chemical potential of the noninteracting system. Due to the smallness of G 
we can assume that the self-energy is also small and solve the equation (1261) by iterations 

/i = /i + ReS(/i, /i). (27) 

Thus, as discussed in Ref. [50] (see also references therein), in the presence of the 
electron-phonon interaction one should take into account the renormalization of the 
noninteracting chemical potential /z to its interacting value fi. In the numerical results 
presented below we fix the value of /i and then in accordance with Eq. (1271) compute the 
value of jl which has to be used as a true chemical potential of the interacting system. 
As one can see from Eq. ( 125|) . for the massless carriers in graphene the role of Ae// is 
to renormalize their velocity: vp vf/{1 + Ag//) and the value of chemical potential 
fi — ^ /i/(l + Ae//) ~ fl. Restoring the imaginary part of the self-energy, F = — ImE 
in (1251) . we find that the quasiparticle width is also renormalized: F F/(l + Xeff)- 
These properties allow us to use previously obtained expressions with the corresponding 
renormalizations. 

Fig. [2^ compares results for the ac conductivity Rea^xi^) as a function of Q to 
800 meV with and without phonons at two temperatures T = 10 K and T = 100 K. For 
T = lOK the dash-dotted (black) curve was calculated using bare bands, a constant 
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CM 

ri.o- 



6jo^0.2eV 

G=0.2 
/;^190 meV 
r=0.5 meV 




Figure 2: (a) The conductivity Keaxxi^) in units 7re^/2/i (universal background) for 
bare band: dash-dotted (black) curve for T = 10 K and short dashed (green) line for 
T = 100 K. The curves with phonons {ujq = 0.2 eV and G = 0.1) are the long dashed 
(red) line for T = 10 K and the solid (blue) line for T = 100 K. In all cases the chemical 
potential /i = 190 meV and the scattering rate F = 0.5 meV. (b) The conductivity with 
phonons same as (a), except to G = 0.2, at various temperatures: T = lOK - the long 
dashed (red) line, T = 50 K - the solid (blue) line, T = 100 K - the dash-dotted (black) 
line and T = 200 K - the short dashed (green) line. 



residual scattering rate F = 0.5 meV and a chemical potential /i = 190 meV. This curve 
is to be compared with the long dashed (red) curve which includes phonons with G = 0.1, 
ujo = 0.2 eV and is at the same temperature. To regularize computation when the self- 
energy ImE(co',/i) given by f|T9l) is zero we use the same value of the residual scattering 
rate F = 0.5 meV. While the main rise in the conductivity to its universal background is 
at 2fi = 380 meV in the bare band case, the electron-phonon renormalization according 
to Eq. fl271) have shifted this rise to 2fl, the renormalized chemical potential which falls 
at considerably lower energy. Note also that at higher energies the bare conductivity 
has risen to its universal value, while the renormalized one falls below, i.e. is slightly 
less than one in the chosen units. Increasing the coupling G reduces the background 
value even more as can be seen in Fig. 2 of Ref. [28]. We find that for G = 5 (which 
is certainly larger than expected for graphene) the conductivity aX Q = 800 meV has 
dropped to a value of 0.87. This reduction is expected on physical grounds. Phonons 
shift optical spectral weight to energies beyond the bare band cutoff on the scale of a 
few phonon energies. This spectral weight must come from lower energy to conserve the 
optical sum. We find the depletion starts for Q > and shows no structure. There 
are no wiggles (excess conductivity associated with the Stokes and anti-Stokes lines) of 
the kind reported in Refs. [26,27]. Our result agrees with an earlier observation made 
in Ref. [28] that these features disappear when the real part ( !20!) of the self-energy is 
included. This is further confirmed in [51]. The reduction in universal conductivity 



On the universal AC optical background in graphene 



11 



found here is important, because it shows that the ac conductivity of graphene does not 
necessarily give metrologically accurate value of the von Klitzing constant h/e"^, although 
it does give a good estimate. A second feature of the long dashed (red) dressed curve as 
compared with dash-dotted (black) curve for bare bands is the region starting at fl = ujq 
and extending to around 2/i (twice the dressed chemical potential). The bare case is 
totally flat and small in this region, but the dressed curve shows additional absorption 
starting with a clear "kink" ai Q = ujq and increasing substantially as Q grows towards 
2/i. These are the Holstein boson assisted absorption processes which proceed with the 
creation of a hole-particle pair and in addition a phonon is present in the final state. 
These processes are well known in ordinary metals and correspond to the incoherent part 
of the electronic spectral density brought about through coupling to the phonons. At the 
same time the coherent Drude like contribution centered at f2 = is reduced by a factor 
1-l-Ae// and contains only a fraction 1/(1-|-Ae//) of the optical spectral weight, while the 
remainder Ae///(1 + Kff) is to be found in the incoherent phonon assisted Holstein side 
band [52]. The short dashed (green) and solid (blue) curve give results at T = 100 K. 
First note that temperature does not significantly change the value of Kea^xi^) at large 
fl and so the ac background value is robust with respect to temperature changes. In the 
region of the main rise around the bare (for the short dashed curve) and renormalized 
(for solid curve) chemical potential, however, temperature does provide smearing. For 
the bare band this smearing is more symmetric about twice the chemical potential than 
it is for the phonon renormalized case, where the signature of a phonon contribution 
remains as extra absorption in the region above the phonon energy uo and below 2/i. 
Comparison with the the long dashed (red) curve at T = 10 K with that at T = 100 K 
shows that, the phonon structures are themselves smeared by temperature and are best 
seen at low temperature. 

This is seen better in Fig. |2]3 where we show the phonon renormalized value of 
Rea^xi^) for four temperatures, namely, T = 10 K, 50 K, 100 K, 200 K. By T = 200 K 
the short dashed (green) curve has become rather smooth and no sharp phonon structure 
remains, but smearing about = 2/i is nevertheless still quite asymmetric and in that 
sense retains information on phonon renormalization. 

There is another way to see phonon renormalization in lieaxxi^) curves that is 
worth commenting on. It was shown in [50] that at T = for |/i| > Q (photon energy) 
and small residual scattering rate F the low frequency conductivity is dominated by 
the intraband transitions and that Keaxxi^) takes on a Drude like profile with a single 
change over the bare band case, namely, 

4r 

with the phonon renormalization Ae// as defined in Eq. fl2T!) and in the lowest 
approximation given by Eq. ( l22l) . Defining the ratio of renormalized to unrenormalized 
conductivity as R{Q, T) we get at zero temperature 
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which starts at value 1 for = and saturates to a value of 1/(1 + \effY at 2V. 
Numerical results for a value of R{VL, T) are shown in Fig. [3] for four temperatures 
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Figure 3: The ratio of dressed (by phonons) to bare conductivity R{Q) defined in the 
text vs fl in the low-energy (to 6 meV) at four temperatures: T = 10 K - the long dashed 
(red) line, T = 100 K - the solid (blue) line, T = 200 K - the dash-dotted (black) hue 
and T = 300 K - the short dashed (green) line. 

T = 10 K, 100 K, 200 K, 300 K. All curves start at value 1 and are close to a saturation 
value of ~ 0.72 by = 6meV with a small decrease in saturation value at higher 
temperature taken to mean that the effective phonon renormalization Ag// has increased 
slightly with temperature. This is a well known effect in simple metals [52]. Of course, in 
an actual experiment one does not have an access to the bare band conductivity and the 
ratio R{Q, T) cannot be formed. However, one can take ratios at different temperatures 
to verify that the electron-phonon renormalization is nearly a constant in temperature 
as found in this work. What is more important is that Fig. [3] shows clearly that the 
intraband Drude type renormalizations of Eq. ( l28l) are significant for the parameters 
used here and should be observable. 

While we have digressed somewhat from our main theme, to describe more general 
effects of phonons on the conductivity, our important result is that, in principle, the 
electron-phonon interaction will reduce the value of the universal background below its 
bare band value. Admittedly for the parameters characteristic of graphene this effect 
is not expected to be large. We note that in Ref. [53] the conclusion was made that in 
general the effects of Fermi liquid corrections on the universal background are small. 

5. AC conductivity in an external magnetic field without phonons 

As we already mentioned in the Introduction, in the papers [21-23,36] we have derived 
and analyzed various expressions for the diagonal and off-diagonal magneto-optical 
conductivity of graphene. The derivation is also based on the analytical calculation 
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of the bubble diagram ([H]), but instead of the GF one has to substitute the GF 
describing Dirac fermions in an external magnetic field. 

Particularly useful are the Eqs. (11) and (12) of Ref. [23]. The first of these 
equations is for the complex diagonal conductivity 

e^vl\eB\{n + 2iT) 



oo 



X 

n=0 



TTCt 



[n^(M„) - n^(M„+i)] + [n^(-M„+i) - n^(-M„)] 



A2 



[npi-Mn 



(30) 



Mn+1 + M„ 



where the energies of the relativistic Landau levels are 



±M„ 



A'^ + 2nvl\eB\/c. (31) 

In Eq. ( !3Ti) the Landau energy scale i^(-B) associated with the magnetic field expressed 
in the units of temperature and energy reads [restoring fi\ 



L\B) 



2eBvl 



^^^fi = 1.7696 X 10-^K24(m/s)5[T], 
ck^ 



(32) 



L2(5) ^ 1314meV25[T], L{B) ^ 36.3meVV5[T], En^±VnL{B), 



where the field B is given in Tesla and vp is the Fermi velocity in graphene given in m/s 
and in the second line it is chosen to be Wi;' ~ 10^ m/s [in what follows we set A = 0]. It 
is assumed in Eq. ( !30|l that all Landau levels have the same width, for the case of the 
Landau levels with different widths see Eq. (9) of [23]. 

To study the limit of a weak magnetic field, we apply the Poisson summation 
formula and express the conductivity fl30l) as follows 



4:71 h 



n k=l n 



(33) 



where f{x) = /*"*''°(x) + /*"*'^^(x). The intraband and interband contributions are given 
by 



intra , 



X) 



npiy/x) — rip{yx + lJ) + rip{—yx + I?) — np^—y/x) 



^ X ^ L? — y/x 



X 



2r 



+ 



2r 



(V^TX^- v^-fi) +(2r)2 (V^TX^-y^ + r]) +(2r)^ 



(34) 



and 

jinter j 



X] 



np{—y/x) — rip{yx + I?) + rip{—yx + I?) — np{y/x) 



\/x + L'^ + y/x 



2T 



+ 



2r 



(V^TX^+v^-fi) +(2r)2 (v/^TX^+ v^ + fi) +(2r)^ 



(35) 
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Although the representation (155]) - (|53]) looks rather complicated, its limiting case 
provide a lot of insight on the behaviour of the conductivity. For example, evaluation 
of Eq. ( l33l) in the high frequency limit, fl ^ L{B), \fi\,T, gives 

Rea^^{il) ~ (36) 



J 41^;^ 



A:=l 



1 



+ 



_(M + l)2 + (2r/fi)2 (n- 1)2 + (2r/fi)2 



COS I vrA; 



''2L^{B) ) 



The first term of Eq. fl5Bl) corresponds to the universal ac background, while the second 
term describes an oscillatory behavior about that value. In Fig. H] a we illustrate this 
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Figure 4: (a) The conductivity Keaxxi^) in units TTe^/2h (universal background) for 
temperature T = 10 K, scattering rate F = 2.5 meV. Long dashed (red) line, the 
chemical potential /x = 30meV and the magnetic field B = 0.1 T, dash-dotted (black) 
line /i = 30meV and B = 0.2 T, solid (blue) line /i = 30meV and B = 0.3 T, short 
dashed (green) line fi = 20meV and B = 0.2 T. (b) The conductivity R.eaxx{^) 
in units ire^ /2h (universal background) for magnetic field B = 0.1 T, temperature 
T = 10 K, and scattering rate F = 1 meV. Long dashed (red) line, the chemical potential 
/i = 50meV, dash-dotted (black) /i = 40meV, solid (blue) fi = 30meV, short dashed 
(green) fi = lOmeV. All thick lines are computed using the full expression (130|1 and the 
thin lines using the first line of Eq. (|39l) . 



analytical result and present the ac conductivity Keaxxi^) in units ire^ /2h as a function 
of VL to 150 meV computed using the formula (1501) . It turns out that even for the smallest 
value of the field B = 0.1 T considered it is sufficient to include 100 terms in the sum. 
We observe indeed that for large Q the conductivity oscillates around the universal 
value. 

In the clean limit, F ^ 0, and for fields L{B) ^ |yu| the main contribution to the 
integrals in Eq. (|33|) comes from the Lorentzians, because they reduce to delta functions, 
and thus we obtain 



axx{^) = ^ <{ [npi-L] - npiL)] 



2TL 



2TL 



[7r[(L-r])2 + (2F)2] Tx[{L + QY + {2Ty 
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+ 



X 



Up 



2n 



Up 



— V 



X COS 



k=l 



dx 



2n 



2T 



Up 



2VL 



+ np 



2n 



2r 



(37) 



(v/^TX2-7i-fi) +(2r)2 (7^TX2 + yi-^]) +(2r)2_ 



(^)]} 



where we kept finite F in the places crucial for the regularization of the conductivity. 
This high field regime L{B) ^ is, in principle, well described by the original 
representation (l30l) . For example, the first term of fl37|) corresponds to the absorption 
line at Q = El — Eq = L{B). In the weak fields, L{B) < the intraband part 
has an additional sharp (at small temperatures) maximum which gives the following 
contribution 



dxf{x) 



dx 



dnp{y/x) ^ dnp{—y/x) 



X 



2r 



2^ 



+ 



2F 



(3^ 



Since the maximum of the integral is at the point x = /i^, one can replace E^ l2^fx by 
the cyclotron frequency uj,. = \eB\vp/{c\iJ,\j^, then the remaining integrals are evaluated 
exactly and we obtain 



Rea,,(fi) = ^Tln (2 cosh 



2F 



+ 



+ 



2h 
'Ah 



[np{-E) - np{E)] 



2F 



+ 



2F 



tanh 



fi + 2/i 



+ tanh 



(L-fi)2 + 4F2 

n - 2fi' 



+ 



2F 



(L + ^])2 + 4F2 
Ei = L(5)<|/i|, 



(39) 



4T 4T 

where similarly to Eqs. ([T]) and fll4p . the impurity scattering rate F was set to zero in 
the interband term. We note that setting 5 = in Eq. fl39|) . we arrive at the Eqs. (13) 
and (14) from Ref. [23] mentioned in Sec. 13.11 Eq. fl39|) illustrates that the universal 
ac background is robust to the perturbations provided by a weak magnetic field. For 
nonzero, but small B such that Ei = L{B) < (or equivalently, for 2/ia;c < the 
first line of (1391) with the same coefficient In [2 cosh(/i/2T)] as in the B = expression 
f|T^ . is much larger than the second line of fl39l) . In this case, the spectral weight 
Wintra givcu by Eq. ( |T6i) of the Drude peak at f2 = from Eq. (IT^ is transferred to the 
cyclotron peak at Q = Uc which has exactly the same field independent spectral weight. 
Further, as the field increases, the dependent on this field spectral weight of the second 



I Accordingly, we have huJc[meY] = 657nieV'^i?[T]//i[meV] 
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line also grows and a second peak a.t Q = Ei develops. As the field grows and becomes 
El = L{B) > one has to use the appropriate in this limit expression fl37j) which has 
only one low-energy peak aX Q = Ei. 

In Fig.Hlb we show the evolution of the Drude peak to the cyclotron peak at = tUc 
observed when Ei < \fi\. One can see that for /i = 50meV the thick long dashed (red) 
is undistinguishable from the thin long dashed (red) computed using the approximate 
expression based on the first line of Eq. fl39l) . For fi = 40meV thick and thin dash- 
dotted (black) lines are already distinguishable, but show that the agreement between 
the full expression (l30l) and the first line of Eq. (l30l) is very good. The agrement remains 
good for lower value of = 30 meV [for solid (blue) lines] . However, when /i is getting 
close to El = 11.5 meV as shown in the short dashed (green) curve for fi = 12meV, 
the approximation based on ( !39l) loses its validity. The position of the first peak on 
the thick short dashed line is shifted from Q = Uc to lower energy and second peak at 
Q = El develops. 

6. Discussion 

We have studied how the universal ac conductivity background in graphene is affected by 
a finite chemical potential, finite temperature, elastic impurity scattering, application 
of a weak magnetic field, as well as the opening of a finite gap which changes the 
quasiparticles into massive relativistic-like fermions. While large changes in the ac 
conductivity are found to arise at low energies, when the photon energy is increased to a 
value much larger than any of the energy scales involved, the corrections to the universal 
conductivity coming from the above mentioned sources are small and the background 
value remains close to TTe'^/{2h). On the other hand, the changes found when the 
electron-phonon interaction is accounted for can be more significant for large values of 
coupling to high energy phonons. These renormalizations increase the band width on 
the magnitude of the phonon energy scale and thus transfer optical weight to states 
beyond the bare band cutoff and consequently reduce the background at intermediate 
energies. Of course, if the actual electron-phonon coupling in graphene is very small 
these effects will also be small. Nevertheless, it is important that the background value 
can be affected by various agencies and need not be exactly Tce'^/{2h). The electron- 
phonon interaction also renormalizes the value of the chemical potential measurable as 
a threshold in the absorption. 

We find that this interaction, while providing phonon assisted sidebands called 
Holstein processes in metal physics, does not introduce wiggles in the background value 
as confirmed in [51]. We demonstrated that in a weak external magnetic field the 
universal value of the ac conductivity does not change while the conductivity itself 
oscillates around this value. Finally, we studied the evolution of the Drude peak in 
graphene with a magnetic field. 
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